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Abstract 



In this talk I will review some of the recent applications of the replica theory to 
glasses. I will firstly describe the basic assumptions and I will show that they can be 
considered as a precise reformulations of old ideas. The relation of this approach with 
the mode-coupling theory will be shortly discussed. 

I will present numerical simulations for binary mixtures. The results of these simula- 
tions point toward the correctness of the replica approach to glasses. I will describe the 
results of off-equilibrium simulations for large systems, in which the aging dynamics is 
studied. 

In this talk I will review some of the theoretical progresses that have been done in our 
understanding of glasses. They have been done mainly by comparing the results obtained for 
soluble models of generalized spin glasses with structural glasses assuming that the phase space 
of the two models are similar. 

In the nutshell many of the ideas I am going to present are not new: they are already in the 
original papers of Gibbs and Di Marzio. However the comparison of glasses and generalized 
spin glasses, introduced in ref. allow us to put these ideas in a sharper form and to test them 
in numerical (and eventually real) experiments. In this talk I will not discuss the theoretical 
basis under which this scenario has been derived (i.e. the mathematical tool needed to derive 
the results stated for the generalized spin glasses) but I will concentrate the attention of the 
physical picture. 

The basic ideas are quite simple [0, H, ||, ||, |^, |^, |]. Let us consider a system of particles 
and let us assume that we can introduce a free energy functional F[p] which depends on the 
density p(x) and on the temperature. We suppose that at sufficiently low temperature this 
functional has many minima (i.e. the number of minima goes to infinity with the number 
(A^) of particles). Exactly at zero temperature these minima coincide with the mimima of the 
potential energy as function of the coordinates of the particles. Let us label then by an index 
a. To each of them we can associate a free energy and a free energy density = Fa/N . 

In this low temperature region we suppose that the total free energy of the system can 
be well approximated by the sum of the contributions to the free energy of each particular 
minimum: 
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When the number of minima is very high, it is convenient to introduce the function 
A/'(/, T, A^) which is the density of minima whose free energy is equal to /. With this no- 
tation we can write the previous formula as 

Z = j dfexp{-pNfmf,T,N). (2) 

In the region where Af is exponentially large we can write 

Ar(/,T,iV)^exp(iVS(/,r)), (3) 

where the function S is called the complexity or the configurational entropy (it is the contri- 
bution to the entropy coming from the existence of an exponentially large number of locally 
stable configurations). 

The relation (H) is valid in the region fm(T) < f < fuiT). The minimum possible value of 
the free energy is given by fm{T). Outside this region we have that J\f{f,T) = 0. It all cases 
known J]{fm{T),T) = 0, and the function S is continuous at fm- 

For large values of we can write 

ffM 

exp{-N(3fs) ^ / dfexp{~N{(3f - S(/,r)). (4) 

We can thus use the saddle point method and approximate the integral with the integrand 
evaluated at its maximum. We find that 

f^fs = mm^f)^Pr-^f\T), (5) 

where 

$(/)^/3/-s(/,r). (6) 

(This formula is quite similar to the well known homologous formula for the free energy ,i.e. 
Pf = miiiEiPE — S{E)), where S{E) is the entropy density as function of the energy density.) 
If we call /* the value of / which minimize $(/). we have two possibilities: 

• The minimum is inside the interval and it can be found as solution of the equation 
P = dTj/df. In this case we have 

^$ = (3f* - S*, S* = S(/*, T). (7) 

The system may stay in one of the many possible minima. The number of minima at 
which is convenient for the system to stay is exp(iVS*) . The entropy of the system is 
thus the sum of the entropy of a typical minimum and of S*, which is the contribution to 
the entropy coming from the exponential large number of microscopical configurations. 

• The minimum is at the extreme value of the range of variability of /. We have that 
/* = fm and ^ = fm- In this case the contribution of the complexity to the free energy 
is zero. The different states who contribute the the free energy have a difference in 
free energy density which is of order A^~^ (a difference in total free energy of order 1). 
Sometimes we indicate the fact that the free energy is dominated by a few different 
minima by say the the replica symmetry is spontaneously broken p!0[| . 



Form this point of view the behaviour of the system will crucially depend on the free energy 



landscape |11], i.d. the function S(/, T), the distance among the minima, the height of the 
barriers among them... Although our final task should be to compute the properties of the free 
energy landscape from the microscopic form of the Hamiltonian, we can tentatively assume 
that the landscape of fragile glasses is similar to that of some soluble long range models in 
presence of quenched disorder |jl|, . 



I cannot discuss here in details the rationale for this hypothesis; it is also clear that it 
cannot be exact and some differences should be present among the predictions of the mean 
field approximation and the real world. Here I will present the scenario coming from mean 
field, stressing some of the predictions that should have a wider range of validity and comparing 
them with numerical simulations for fragile glasses. 

We can distinguish a few temperature regions. 

• For T > Tf the only minimum of the free energy functional is given by a constant density. 
The system is obviously in the fluid phase. 

• For Tf > T > To there is an exponentially large number of minima with a non-constant 



density p{x) |T2| . There are values of the free energy density such that the complexity 
E is different from zero, however the contribution coming from these minima is higher 
that the one coming from the liquid solution with constant p{x). 

The most interesting situation happens in the region where Td > T > Tk- In this region 
the free energy is still given the fluid solution with constant p and at the same time the 
free energy is also given by the sum over the non trivial minima [|^, 

Although the thermodynamics is still given by the usual expressions of the liquid phase 
and final free energy is analytic at T^, below this temperature the liquid phase correspond 
to a system which at each given moment may stay in one of the exponentially large 
number of minima. It is extremely surprising that the free energy of the liquid can be 
written in this region as the sum of the contribution of the minima, according to formula 

i)- 

The time to jump from one minimum to an other minimum is quite large: it is an 
activated process which is controlled by the height of the barriers which separate the 
different minima. The correlation time will become very large below and for this 
region is called the dynamical transition point. The correlation time (which should 
be proportional to the viscosity) should diverge at T^- The precise form of the this 
divergence is not well understood. It is natural to suppose that we should get divergence 
of the form exp(y4/(T — T^Y) for an appropriate value of z/, whose reliable analytic 
computation is lacking [|T], 

The equilibrium complexity is different from zero (and it is a number of order 1) when 
the temperature is equal to To and it decreases when the temperature decreases and it 
vanishes linearly at T = Tk- At this temperature (the so called Kauzmann temperature) 
the entropy of a single minimum becomes to the total entropy and the contribution of 
the complexity to the total entropy vanishes. 

In the region where T < Tk the free energy is dominated by the contribution of a few 
minima of the free energy having the lowest possible value. Here the free energy is no 
more the analytic continuation of the free energy in the fluid phase. A phase transition 
is present at Tk and the specific heat is discontinuous here. 



o 



Figure 1: The qualitative dependence of the free energy as function of the configuration space 
in the region relevant for the dynamical transition, i.e. for T < T^. 



The free energy landscape is quit usual; we can try do present the following pictorial 



interpretation, which is a rough simplification |13|. At temperature higher than To the system 



stays in a region of phase space which is quite flat and correspond of a minimum of the total 
free energy. On the contrary below the phase space is similar to the one shown pictorial 
in fig. 1. The region of maxima and minima is separated by the region without barriers by 
a large nearly flat region. The minima in the region at the left are still present also when 
Tf > T > Td, but they do not correspond to a global minimum. 

At temperatures higher than To the system at thermal equilibrium states in the right 
region. When the temperature reaches the system arrives in the fiat region. Here the 
potential is fiat and this causes a more or less conventional Van Hove critical slowing down 
which is well described by the well known mode coupling theory (which is exact in the 
mean field approximation). The mode coupling theory describes the critical slowing down 
which happens near [^. 

In the mean field approximation the height of the barriers separating the different minima 
is infinite and the temperature To is sharply defined as the point where the correlation time 
diverge. In the real world activated process (which are neglected in the mean field approxima- 
tion and consequently in the mode coupling theory) have the effect of producing a finite (but 
large) correlation time also at T^. The precise meaning of the dynamical temperature beyond 
mean field approximation is discussed in details in 

When the temperature is smaller that we must be more precise in describing the dynam- 
ics of the system. Let us start from a very large system (of particles) at high temperature 
and let us gradually cool it. We find that it should go at equilibrium in the region with many 
minima. However coming from high free energy (from the right) it cannot enter in the region 
where are many maxima and minima if we wait a finite amount of time (the time to crosses 
the barriers diverges as exp{AN). If we do not wait an exponentially large amount of time the 
system remains confined in the fiat region. In this case the so called dynamical energy, 

Ed = lim lim E{t,N), (8) 

is higher that the equilibrium free energy. The situation is described in fig. 2. 

Below Td the system is trapped in metastable states when cooled. The time needed to 
escape from these states diverges when A^ goes to infinity in the mean field approach where 
activated processes are forbidden. Of course the difference of the static and dynamic energy is 




Figure 2: The qualitative behaviour of the equihbrium and of the dynamical energy as function 
of the temperature. 



an artifact of the mean field approximation if we take literarily the limit t — oo in the previous 
equation because as matter of fact there are no metastable states with strictly infinite mean 
life. However it correctly describe the situation on laboratory times, where metastable states 
are observed. 

In real systems, beyond the mean field approximation, the height of the barriers is finite 
also below To and the mean life of the metastable states is finite, albeit very large. Some 
mechanisms have been described who imply a divergence of the correlation time in real sys- 
tems at the Kauzmann temperature 0, |^, |16|. The conventional glass temperature, i.e. the 
temperature at which the microscopic correlation time becomes macroscopic (e.g. of order of 
the minute) is between the two temperatures (i.e. Td > To > Tk). 

It should be clear that in this framework the dynamical temperature Td is not so well 
defined and it correspond to a crossover region below it the dynamics is dominated by activated 
processes 0. 

In the mean field approximation there are very interesting phenomena that happen below 
Td when the system is cooled from the high temperature phase. These phenomena are related 
to the fact that the system does not really go to an equilibrium configuration but wanders in 
the phase space never reaching equilibrium. The phenomena are the following: 

• The energy approaches equilibrium slowly when the system is cooled from an high energy 
configuration. In other words for large times we have 

E{t,T) = Ed{T) + B{T)r^^'^\ (9) 

where the exponent A(T) does not vanish linearly at zero temperature as it should hap- 
pens for an activated process. 

• Aging is present, i.e. the correlation functions and the response functions in the region 
of large time do depend on the story of the system |l^ . 



In the region where aging is present the hypothesis at the basis of the fluctuation dissi- 



pation theorem are no more valid. New generalized relations are satisfied ||20| , pT| , which 
replace the fluctuation dissipation theorem. 

In the region where the diffusion constant is zero, a new phenomenon, logarithmic diffu- 
sion, is present. 



All these phenomena are well known also in the case of spinodal decomposition but have a 
more general validity. 

Here I will illustrate this phenomena in the case of a monatomic fluid (I will consider the 
case of a binary mixture in order to avoid crystalization). The particle interacts with a soft 
potential. Half of the particles are of type A, half of type B and the interaction among 
the particles is given by the Hamiltonian: 



a{i) + a{ky 
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i<k \ -^^'l / 

where the radius (a) depends on the type of particles. 

The model is well studied in the past ^ |2^, |2^. It is known that the choice 
o'b/o'a = 1-2 strongly inhibits crystallisation and the system goes into a glassy phase when it 
is cooled. Using the same conventions of the previous investigators we consider particles of 
average radius 1, i.e. (cr^ + 2((j^ + ctb)^ + ct|)/4 = 1. The system goes into the glass phase at 
r ^ 1.45, where we follow the standard notation: T = (3^/^. 



For reason of space I will concentrate on the properties of diffusion |]26|. Let me firstly 
discuss the equilibrium situation and later I will discuss the non-equilibrium case. 

A typical quantity that can be studied in a glass is the diffusion. We can define the average 
distance squared done by a particle in time t as 

m = {Mt)-xmf) (11) 

where the label i may take values from 1 to A^. The diffusion constant D is defined as 

A(t) ^ Dt (12) 

for large t. 

In the following we will assume (and it is a rather good approximation) that in the glass 
phase A(t) goes to a constant (A*) at large t, so that the diffusion constant D is zero. 

We can also introduce the response to an external force. Let us consider the following time 
dependent Hamiltonian 

H = Ho + eit)F ■ Xi, (13) 

where Hq is the original unperturbed Hamiltonian and F is a vector of length where d is 
the dimension of the space (usually 3). Both here and in eq. (0) the value of index i (i.e the 
particle we have chosen) is arbitrary and the result does not depend on i. 
We consider the case where 



e(t) = for t < 0, e{t) = e for t > 0, (14) 
e being a small number. The response function is defined as 

{F ■ x,{t)) = eRit) + {F ■ x,{t))\e=o- (15) 
The usual fluctuation dissipation theorem tell us that 

(3m = m, (16) 

In the large t limit we obtain the celebrated Einstein relation among diffusion and viscosity. 

Let us consider an off-equilibrium situation. Here the results depend on the history of the 
system. We a very simple history. The system is in the liquid phase for t < 0, at t = it is 




cooled very fast at the final temperature in the glassy phase. Here time translational invariance 
is no more present. We generalize the previous definitions introducing two times quantities: 

A{t,t^) = {ixiit + t^)-Xi{t^)y). (17) 

The response function is defined by 



eR{t, t^) = {x,{t + ■F)-{F-x^{t + t^)) 



e=0 



where the statistical expectation values are computed with the time dependent Hamiltonian 
(O), where e{t) given by 



e{t) = for t < t^, e{t) = e hi t > (19) 

We obviously have that 

lim A{t,t^) = A{t) (20) 

and in this way we recover the equilibrium limit. 

The equilibrium limit is relevant when » t. The aging regime is present when both t 
and tw are both large. In the simpler case (i.e. naive aging) the interesting behaviour happens 
when t = 0(tw) (i.e. the ratio t/tw is finite) and is large. In this regime we find that 

A(t,t^)^ a(-^). (21) 

In other words if we measure something on a time scale t comparable to t^^ we always see a 
non-equilibrium behaviour. 

The generalized fluctuation dissipation theorem tells us that, if we plot at a given A 
versus i?, the plot has a limit when t^, goes to infinity (as an effect of aging) and we can define 
a function X(A) as 

(9/? 

^ = /5^(A). (22) 

The function X can be constructed for any pairs of observable and its value should not 
depend on the observables. More precisely, if we take two other quantities satisfying the 
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Figure 4: The quantity m = dR/dj3A as = 2048 as function of the temperature. The 
straight hne is the prediction of the approximation m(T) = T/T^. 



fluctuation dissipation theorem, we should have that at equal values of s the function X takes 
the same value of X{A). The fluctuation dissipation theorem is violated in the same way for 
all the observables at the same time. 

It has been conjectured that this dynamically defined function X is related to an similar 
function which can be defined by measurements done only at equilibrium p|. A precise presen- 
tation of this conjecture would take too much time. A numerical verification of its correctness 



in short range spin glasses can be found in [21 



The simplest non trivial form of the function X(A) is (in analogy to what happens in 
generalized spin glasses) 



X(A) = 1 for A < A* 
X(A) = m(T) for A > A* 



(23) 



where m(T) is approximately a linear function of the temperature taking the value one at 
T = Td. 

In other words if we plot the response R as function of (3 A we should find piece-linear graph 
with slope 1 for at small A and slope m for large A. The result of numerical simulations in 
the glassy phase are shown in fig. 3 (for details see ref. 



Naive aging would imply that in the large time region (both times large) A(t,t^) is a 
function of only the ration t/t^. If we do some reasonable assumption naive aging implies that 
/S.{t,tw) behaves as A\n.{t/tu]) + B when t is large [p^ . 

This phenomenon is well known in Ising spin systems. It amount to say that if we cool 
at zero temperature a spin system starting from an high energy configuration, the average 
numer of times a spin is reversed increases as the logarithm of the time. This off-equilibrium 
logarithmic diffusion seems to be present in glasses. For example numerical simulations of a 
binary mixture are very well in agreement with this prediction as can be seen in fig. 5. 

The results of simulations are in very good agreement with the theoretical expectations 
based on our knowledge extracted from the mean field theory for generalized spin glass models. 
The approximation m(T) = T/Tq seems to work with an embarrassing precision. We can 
conclude that the ideas developed for generalized spin glasses have a much wider range of 
application than the models from which they have been extracted. It likely that they reflect 
quite general properties of the phase space and therefore they can be applied in cases which are 
very different from the original ones. In a recent work 0] some thermodynamic predictions have 
been obtained for the behaviour of glassy systems, like the present ones, under the assumption 
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that they obey the laws derived for generahzed spin glasses. The present work confirms that 
assumption. 

The most urgent theoretical task would be now to develop an analytic theory for glasses 
in the low temperature region from which one could compute the function m(T). This goal 



should not be out of reach: a first step in this direction can be found in p9 
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